Exchange interactions and high-energy spin states in Mni 2 -acetate 
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We perform inelastic neutron scattering measurements on the molecular nanomagnet Mni2-acetate 
to measure the excitation spectrum up to 45 meV (500 K). We isolate magnetic excitations in two 
groups at 5-6.5 meV (60-75 K) and 8-10.5 meV (95-120 K), with higher levels appearing only at 27 
meV (310 K) and 31 meV (360 K). From a detailed characterization of the transition peaks we show 
that all of the low-energy modes appear to be separate S = 9 excitations above the S = 10 ground 
state, with the peak at 27 meV (310 K) corresponding to the first S — 11 excitation. We consider a 
general model for the four exchange interaction parameters of the molecule. The static susceptibility 
is computed by high-temperature series expansion and the energy spectrum, matrix elements and 
ground-state spin configuration by exact diagonalization. The theoretical results are matched with 
experimental observation by inclusion of cluster anisotropy parameters, revealing strong constraints 
on possible parameter sets. We conclude that only a model with dominant exchange couplings 
Ji ~ Ji ~ 5.5 meV (65 K) and small couplings J3 ~ J4 ~ 0.6 meV (7 K) is consistent with the 
experimental data. 



PACS numbers: 75.30.Et, 75.50.Xx, 78.70.Nx 

I. INTRODUCTION 

The topic of molecular magnets 1 ' 2 has emerged in the 
last decade as one of the major interdisciplinary fields of 
research in the materials science community. Molecular 
magnets, also called "spin clusters", are crystalline ma- 
terials composed of magnetic centers, mostly transition- 
metal ions, which have strong mutual interactions within 
each molecule. Each spin cluster in the lattice is magneti- 
cally well isolated from its neighbors due to the presence 
of surrounding ligands. This magnetic shielding allows 
the study of the individual behavior of nanometer-scale 
magnetic systems. 

Chemists and physicists have combined efforts in an 
increasingly coordinated approach to develop and design 
new magnetic molecules which can show purely quantum 
properties at the macroscopic level. Such phenomena in- 
clude very slow relaxation of the magnetization below 
the "blocking" temperature Tb. 3,4 and quantum tunnel- 
ing of the magnetization vector through an energy bar- 
rier between "spin up" and "spin down" configurations. 5 
The barrier A is governed, to leading order, by two pa- 
rameters, the ground-state spin value S and the axial 
anisotropy D < of the cluster, which yield A = —DS 2 . 
The magnetic properties of the molecules are deter- 
mined primarily by the exchange interactions between 
individual ionic spins and by their magneto-crystalline 
anisotropy. Competition among the different exchange 
couplings can lead to a wealth of situations ranging from 
a nonmagnetic ground state (S = 0) to a variety of large- 




FIG. 1: (Color online) Structure of Mni2-acetate viewed along 
the c axis. Large solid spheres represent Mn 3+ ions (outer 
ring) and Mn 4+ ions (inner core). All Mn atoms have a dis- 
torted octahedral coordination geometry. 



spin ground states (S > 10). A precise knowledge of the 
values of the exchange interactions, and their dependence 
on the geometry of the molecular cluster, is therefore crit- 
ical to understand the ground state and, more generally, 
the possibilities for new synthetic procedures providing 
better materials with larger energy barriers. 

Inelastic neutron scattering (INS) is a powerful 
experimental probe of magnetic excitations and ex- 
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change or anisotropy parameters in molecular magnetic 
clusters. 6-15 Here we investigate the high-energy mag- 
netic states in the best-characterized molecular magnetic 
system to date, Mni 2 -acetate, using INS in the energy 
range between 1 and 100 meV (10-1000 K) and with 
no external magnetic field. Comparison between our ex- 
perimental results and numerical calculations enables us 
to estimate the exchange interaction parameters in the 
Mni2-acetate cluster. 

The manuscript is organized as follows. In Sec. II we 
review the structure and properties of Mn 12 -acetate, ex- 
perimental results, and attempts which have been made 
to understand the exchange interactions. In Sec. Ill we 
present the results of INS measurements performed on 
three different spectrometers, and use these to charac- 
terize the magnetic excitations. Section IV contains a 
theoretical analysis of the conventional exchange model 
for Mni2-acetate by high-temperature series expansion 
techniques for the susceptibility and exact diagonaliza- 
tion to obtain the ground state and low-lying excited 
states, both refined by comparison with the existing data. 
A summary and conclusions are presented in Sec. V. 



II. STRUCTURE AND PROPERTIES OF 
Mni2- ACETATE 

A. Structural Information 

Mni2-acetate 4> 16-18 is a mixed- valence (Mn 3+ /Mn 4+ ) 
compound where the magnetic ions are arranged in two 
groups: a central core composed of a tetrahedron of 4 
Mn 4+ ions (S = 3/2) and an external ring, or crown, of 
8 Mn 3+ ions (S = 2). Figure 1 shows the Mni 2 -acetate 
cluster viewed along the c axis. The point group of the 
Mni2 molecule in the crystal structure is S4. To sim- 
plify the analysis we make the additional assumption of 
fourfold rotation and reflection symmetry of the individ- 
ual molecular clusters, and return later to a discussion 
of this approximation. Each cluster is only very weakly 
coupled to its neighbors, which are separated from each 
other by molecules of water and acetic acid, such that 
no long-range magnetic order has been found for tem- 
peratures as low as the mK range. Consequently, most 
of the experimental work performed at temperatures ex- 
ceeding 1 K may be interpreted in terms of the proper- 
ties of a single molecule. Neighboring Mn ions within 
a cluster are coupled in an intricate pattern by differ- 
ent types of /x-oxo bridge and by acetate bridges, as a 
result of which both antiferromagnetic (AFM) and ferro- 
magnetic (FM) exchange interactions may be present in 
the system. A schematic representation of the exchange 
couplings is shown in Fig. 2. Within the approximation 
of fourfold cluster symmetry there are three incquivalcnt 
Mn sites with four different exchange couplings between 
neighboring Mn ions: J\ (involving two /i-oxo bridges) 
and J2 = Jia ~ -hb (one /z-oxo bridge) between Mn 3+ 
and Mn 4+ ions, J3 between Mn 4+ ions inside the core 




FIG. 2: Schematic representation of the magnetic exchange 
paths in Mni2-acetate as first proposed by Sessoli et al. 16 Ji 
couples Mn 3+ and Mn 4+ ions (d — 2.771 A), J2 a ,b correponds 
to Mn 3+ /Mn 4+ (d = 3.45 A), J 3 to Mn 4+ /Mn 4 + (d = 2.817 
A) and Ji a ,b to nearest-neighbor Mn 3+ ions (d — 3.33 and 
3.41 A). 

tetrahedron (two /i-oxo bridges) and J4 = J^ a « J^h be- 
tween Mn 3+ ions around the external ring (one /i-oxo 
bridge and two carboxylate groups). Inspection of the 
Mn-Mn distances and Mn-O-Mn angles presented in Ta- 
ble I suggests that the approximations J^ a ~ -hb and 
J 4a rs J Ab are eminently reasonable. 



B. Ground State and Anisotropy 

The magnetic ground state has a total spin S — 10, as 
first suggested by Sessoli et al. le ' 17 and later confirmed by 
many experimental studies, including magnetization, 19 
nuclear magnetic resonance (NMR), 20 electron param- 
agnetic resonance (EPR), 21,22 and INS 7 measurements. 
The S = 10 ground state may be rationalized by con- 
sidering a fcrrimagnetic arrangement of 8 parallel spins 
S = 2 on the crown Mn 3+ ions oriented antiparallel to 
the 4 aligned S — 3/2 spins on the core Mn 4+ ions. 

While all of the Mn atoms have a distorted octahe- 
dral coordination geometry, the distortion is significantly 
stronger for the crown Mn 3+ ions. The strongly nega- 



TABLE I: Mn-Mn distances and /i-oxo bridge angles for the 
different exchange couplings in Mni2-acetate. Mn pairs cou- 
pled by one /i-oxo bridge have greater separation (d > 3.3 A) 
and a higher angle than those coupled by two /i-oxo bridges 
(d < 3.3 A). 
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tive axial anisotropy parameter D of the cluster, which 
is due primarily to the Jahn- Teller distortion (axial octa- 
hedral elongation) around the Mn 3+ ions, lifts the de- 
generacy of the S = 10 ground-state multiplet, with 
the M s = ±10 states lying lowest. The M s = ±10 
and Ms = states are separated by an energy barrier 
A = —DS% of approximately 65 K (5.5 meV), which 
blocks the thermal relaxation in a typical magnetization 
relaxation experiment. 17,23-25 Surprisingly, the energy 
barrier A in Mni2-acetate remains the highest known, 
despite intense efforts to synthesize new magnetic clus- 
ters. 

Higher-order (rhombic and quartic) terms in the 
anisotropy are responsible for a mixing of wave functions 
between pairs of states of equal \Ms\ on either side of 
the barrier, opening a channel for tunneling through the 
barrier. This effect has been demonstrated and studied 
extensively in recent years by a variety of techniques, in- 
cluding low-temperature magnetometry, 19 ' 26-28 NMR, 20 
and specific heat. 23,29 An excellent determination of the 
quartic terms has been achieved using high-field EPR 21 
and INS. 7 The recent discovery 30-32 of disorder in the 
acetic-acid and water molecules, which further lowers the 
site symmetry of some of the Mni2 clusters present in 
the lattice, has shed additional light on a long-standing 
controversy concerning the physical origin of the experi- 
mentally observed tunneling. 33-37 

With the primary aim of understanding the mag- 
netic excitation spectrum, in this study we will treat 
the anisotropy contributions of single Mn ions using 
only phenomenological axial anisotropy parameters D s 
of the cluster spin S (D above denotes the ground- 
state anisotropy parameter, denoted f° in Sec. IIIF), 
and neglect all anisotropy terms of higher order than 
quadratic in single-ion or cluster spins. Because single- 
ion anisotropies are some two orders of magnitude 
smaller than the exchange couplings, and higher-order 
anisotropies smaller still, this treatment is expected to be 
acceptable. Nevertheless, approximations at the level of 
single-ion anisotropies do limit the accuracy with which 
experimental results may be fitted in a theoretical anal- 
ysis. 



C. Magnetic Excitations 

While the low-temperature properties of Mni2-acetate 
are well described by an isolated S = 10 ground state, to 
date very little is known about the higher excited states 
whose energies are determined directly by the exchange 
interactions (Ji, J 2 , J3, J4). The first experimental infor- 
mation was provided by magnetization measurements at 
very high fields using an explosive compression technique 
which can access magnetic fields in excess of H = 900 
T. 38 ' 39 A series of peaks is observed between 300 T to 
600 T in the field derivative dM/dH, suggesting impor- 
tant steps in the magnetization curve which may be in- 
terpreted in terms of crossovers from S = 10 to S = 11, 



12, . . . , up to the maximum possible spin value of S = 22. 
Further information on the interaction parameters may 
be obtained from fits to the high-temperature part of 
the zero-field susceptibility. However, because the Curie- 
Weiss regime is not established even at 300 K, 16 and 
significant loss of solvent occurs at and above 308 K, it is 
difficult to draw a reliable conclusion from such fits with- 
out a theoretical treatment of the susceptibility which is 
valid also for lower temperatures. Such an analysis has 
not been performed previously, and will be presented in 
Sec. IV below. 

Based on susceptibility and magnetization data alone, 
the unambiguous determination of exchange constants of 
similar magnitude in such ferrimagnetic systems is im- 
possible. INS is a unique spectroscopic technique which 
can gather information over a wide energy range, and is 
thus capable of identifying a number of energy levels suf- 
ficient to evaluate the exchange interactions in complex 
systems such as Mni 2 -acetate. 

Hcnnion et al. , using an approximately 30% deuter- 
ated Mni 2 -acetate sample, showed clearly for the first 
time that important higher-energy transitions were ob- 
servable above the transitions within the ground-state 
multiplet. These authors reported a series of inelastic 
peaks between 4 meV and 10 meV (45-110 K) which 
were ascribed to transitions between the S — 10 ground 
state and S = 9 excited states. In this study we confirm 
the INS peaks observed in Ref. 6, providing more detailed 
information concerning these, and extend our survey to 
higher energies. Comparison between the experimental 
data and the spectra obtained from exact numerical di- 
agonalization of a four-parameter model for the magnetic 
cluster allows us to assign the observed energy levels and 
to determine a consistent set of values for the exchange 
interactions in Mni2-acetate. 



D. Estimates of Exchange Interactions 

The first theoretical attempts to determine the ex- 
change interactions considered a simplified model 16,40-42 
predicated on the assumption of a single AFM exchange 
interaction 3\ much larger than any of the other param- 
eters in the system. In this situation the 4 Mn 3+ -Mn 4+ 
pairs coupled by J\ are locked in their S = 1/2 ground 
state, and the size of the spin space is considerably re- 
duced, namely to 4 spins 5 = 1/2 and 4 spins S — 2. Al- 
though this model provided a qualitative account of the 
available susceptibility data, it was far from providing a 
unique answer. Further, from qualitative considerations 
of exchange interactions, and particularly the proximity 
of the Ji ^-oxo bridging angle to 90°, even this simplify- 
ing assumption lacks justification. In short, a fully gen- 
eral starting point is required, and the absence of detailed 
information concerning the location of the magnetic exci- 
tations precludes any unambiguous determination of the 
exchange couplings. 

Several attempts have been made to evaluate the ex- 
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change interactions by applying numerical techniques, 
specifically exact diagonalization > and different vari- 
ants of density functional theory (DFT) calculations, 45 ' 46 
to a model of the type discussed above. The results 
of these studies were then compared with a selectively 
chosen subset of the available data, taken from high- 
field magnetization, 38,39 a.c. susceptibility, 23,47 high-field 
EPR, 48 and INS 6 measurements, to argue in support of 
a particular parameter set. Raghu et aZ. 43 considered 
the parameters proposed in Ref. 16, and demonstrated 
that these do not in fact yield an S = 10 ground state. 
Variants on this set were tested which did satisfy this 
condition, but we will show below that the very large 
values of Ji [215 K (18.5 meV)] and of J 2 and J 3 [85 K 
(7.3 meV)] are inconsistent with the measured suscepti- 
bility. Regnault et at, 44 focusing on the S = 10 condi- 
tion, obtained J x = 119 K (10.2 meV), J 2 = 118 K (10.2 
meV), J 3 = -8 K (-0.69 meV), and J 4 = 23 K (1.98 
meV), where the magnitudes of the parameters were de- 
termined by fitting the magnetization data. Park et aZ. 46 
obtained a not dissimilar result in an ab initio DFT cal- 
culation, specifically Ji = 115 K (9.91 meV), J 2 = 84 K 
(7.24 meV), J 3 = -4 K (0.34 meV), and J 4 = 17 K (1.47 
meV). In both of these cases, J\ and J 2 are rather close 
in value, in sharp contrast to the assumption underpin- 
ning the simplified, dimerized model. While J\ is much 
smaller than previously suggested, 16 in combination with 
the larger J 2 value these sets also overestimate the high- 
temperature susceptibility quite considerably. The much 
weaker J 3 and J4 interactions are then difficult to deter- 
mine systematically, although both sets suggest that J 3 
may be weakly FM. Finally, Boukhvalov et al. 45 obtained 
a smaller set of values, J\ = 47 K (4.1 meV), J 2 = 26 K 
(2.2 meV), two different core-spin coupling terms J 3 = 30 
K (2.6 meV) and J3 = 10 K (0.86 meV), and J 4 = 7 K 
(0.60 meV). However, from our computations these last 
values appear to give neither a ground state with S = 10 
nor a suitable reproduction of the magnetization. 



III. EXPERIMENT 



A. Experimental Details 



The INS experiments were performed on three differ- 
ent instruments located at three different sources: 1) 
the time-of-flight (TOF) spectrometer IN4 at the Insti- 
tut Laue-Langevin (Grenoble, France) using wavelengths 
A = 1.1 A and 2.2 A, 2) the TOF spectrometer FOCUS at 
the Paul Scherrer Institute (Villigen, Switzerland) using 
the wavelength A = 3.1 A, and 3) the MARI spectrome- 
ter at the pulsed neutron source ISIS (Didcot, UK) using 
neutrons with incident energies between 12 meV and 50 
meV (in this section we quote energies only in meV). IN4 
and MARI are best suited for high energy transfers and 
a wide Q-range, while FOCUS is a cold-neutron spec- 
trometer designed for energy transfers below 5-10 meV. 
Data were collected at several temperatures between 1.5 
K and 100 K, and corrected for the background and de- 
tector efficiency by means of a vanadium reference and 
empty-cell measurements for each incident energy. 



On IN4 and FOCUS, we used a fresh, 3.8 g polycrys- 
tallinc powder sample of fully deuterated Mni 2 -acetate 
placed under helium in a flat, rectangular aluminum slab 
of 3 mm thickness. Full deuteration was achieved by 
systematic substitution of hydrogenated precursors with 
deuterated versions under argon. On MARI, the sample, 
wrapped in aluminum foil, was mounted in a hollow cylin- 
der of diameter 42 mm and height 62 mm. The sample 
thickness was approximately 3 mm. 



It is thus clear that no consensus has yet emerged on 
the exchange coupling values in Mni 2 -acetate. The ob- 
jective of this study is to determine these exchange in- 
teractions based on the present INS results for the pre- 
viously unexplored energy range up to 45 meV (500 K), 
in combination with information obtained from suscep- 
tibility and magnetization data. We will find that it is 
possible to determine an effectively unique set of interac- 
tion parameters capable of explaining consistently all of 
the available information, i.e. the S = 10 ground state, 
the high-temperature susceptibility, the approximate lo- 
cation of the observed magnetization steps, and the mag- 
netic excitations measured by INS. 



On IN4, the 3 He detector banks cover the angular 
range 29 = 13-120°, giving access to momentum trans- 
fers 0.7 < Q < 4 A" 1 at A = 2.2 A and 2 < Q < 9 A" 1 at 
A = 1.1 A. The resolution obtained from a metallic vana- 
dium reference, given as the full-width at half-maximum 
peak height of the clastic line, was T = 0.9 meV at 2.2 
A and T = 3.8 meV at 1.1 A. On FOCUS, the detector 
banks cover the angular range 29 — 10-130°, giving access 
to momentum transfers 0.7 < Q < 3.4 A -1 at A = 3.1 
A. FOCUS is a time- and space-focusing TOF spectrom- 
eter which allows inelastic focusing. At an energy loss 
of approximately 5 meV, the resolution, obtained from a 
metallic vanadium reference, was T = 0.4 meV at 3.1 A. 
On MARI, the detectors cover the angular range 29 = 12- 
135° and momentum transfers up to Q w 8-9 A -1 for the 
highest energy transfer used (50 meV). The elastic res- 
olution ranged from 0.15-0.2 meV at an incident energy 
Ei = \2 meV to 1 meV at E, = 50 meV. 
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B. Magnetic Neutron Cross-Section 

The differential magnetic neutron cross-section for a 
transition |* m ) -> |*„) is 49 
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In this equation N is the number of Mni2 molecules in 
the sample, k and k! are the wavenumbers of the incom- 
ing and scattered neutrons, Q is the scattering vector, 
cxp(— 2W(Q)) is the Debye- Waller factor, hw is the neu- 
tron energy, |>t m ) are the cluster wave functions with 
energies E m , F(Q) is the magnetic form factor, B4 is the 
space vector of the ith Mn ion in the cluster, a and (3 
represent the spatial coordinates x, y, and z, e and m e 
are respectively the charge and mass of the electron, c 
is the speed of light, and 7 = —1.91 is the gyromagnetic 
constant of the neutron. INS selection rules impose that 
the matrix elements in Eq. (1) arc nonzero only when 
AS = S - S' = 0, ±1 and AM S = M s — M' s = 0, ±1 
where the initial and final states are defined respectively 
by (S,M S ) and (S',M' S ). 



C. Energy-dependence 

Figure 3 shows inelastic spectra obtained on MARI at 
T = 8 K at three different incident energies (Ei — 50, 17, 
and 12 meV). Several transitions are observed between 
1.2 meV and 32 meV. At E t = 50 meV [Fig. 3(a)], four 
sets of transitions are visible above a broad background 
produced by incoherent scattering. These sets are labeled 
(I) to (IV), and their peak energies are 5.3, 8.5, 27, and 
31 meV, respectively. The positions of (I) and (II) are 
consistent with the observations made by Hennion et al. 6 

Figure 3(b), corresponding to Ei = 17 meV, shows 
that peaks I and II are composed of several components. 
These are labeled (la), (lb), (Ha), (lib), and (He), and 
from the center lines of Gaussian fitting curves appear 
respectively at energies of 5.15, 6.04, 7.90, 8.64, and 9.73 
meV. While a double-peak structure is clearly established 
for transitions (la) and (lb), it is less obvious for the 
higher-energy transitions. The peaks were fitted using 
Gaussian lineshapes by assuming a constant line width 
T « 0.75 meV for all transitions and a linearly decreas- 
ing slope to account for the background. The width is 
determined from transitions (la) and (lb) and its value 
is maintained for the other transitions. This is again in 
qualitative agreement with Fig. 8 of Ref. 6. 

Figure 3(c) confirms clearly that there is an energy 
splitting between transitions (la) and (lb) but transi- 
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FIG. 3: INS spectra obtained at T = 8 K on MARI at (a) 
Ei = 50 meV, (b) E z = 17 meV, and (c) E z = 12 meV. The 
Q-range is restricted to 1 < Q < 2 A -1 in panels (a) and 
(b) and to 0.25 < Q < 1.25 A -1 for panel (c). Transitions 
are labeled as discussed in the text. In panel (c), (GS1) is a 
transition within the S = 10 ground-state multiplet and (Ph) 
denotes a low-energy phonon mode. 



tions (Ha) and (lib) remain (barely) unresolved. At this 
value of Ei the width of these transitions is between 0.45 
and 0.55 meV. Based on its Q- and T-dependence (not 
shown), the broad peak around 2.65 meV may be at- 
tributed to phonon excitations. In addition, with this res- 
olution it is now possible to observe the transition peak 
within the ground state, labeled (GS1), at fuv w 1.25 
meV, in complete agreement with the results of Mire- 
beau et al. 7 This peak is significantly sharper (r w 0.35 
meV) than the higher-energy transitions, suggesting an 
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FIG. 4: (a) Inelastic spectra at T = 1.5 K, 10 K, and 20 K 
obtained on FOCUS at an incident wavelength A — 3.1 A. 
At 1.5 K, transitions (la) and (lb) are observed at ftto = 5.3 
meV and hu> = 6.2 meV, respectively. The data taken at 
1.5 K are reproduced as a solid line and shifted upwards to 
facilitate comparison with the 20 K data.(b) Inelastic spectra 
at T = 1.5 K, 20 K, and 50 K obtained on IN4 at A = 2.2 A. 



intrinsic line-broadening of transitions (la) to (lib). 



D. T-dependence 
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FIG. 5: Q-dependence at 1.5 K of peaks (la) and (lb) ob- 
tained from the inelastic spectra shown in Fig. 4(a). Data 
from both FOCUS and IN4 are shown after normalization of 
the maximum intensity at Q ~ 1.45 A -1 . 

new magnetic transitions visible by INS and permits an 
estimate of the anisotropy splitting in the excited states 
(see Sec. Ill F). 

Figure 4(b) presents inelastic spectra obtained on IN4 
over a range of temperatures between 1.5 K and 50 K, 
where the Q-range is limited to 1 < Q < 2 A -1 . At 1.5 K 
one finds the transitions discussed above, which weaken 
progressively as the temperature increases, and finally 
merge into a single, broad feature in which both magnetic 
and phonon scattering are present. A general shift of 
intensity to lower energies is observed with increasing 
temperature for both peaks groups I and II. For the low- 
temperature data, by keeping the line width fixed to the 
value obtained by fitting transitions (la) and (lb), T = 
0.77 meV, the best fit to the large band at higher energies 
is obtained by considering four peaks, (Ha) to (lid), with 
energies of 8.1, 8.8, 9.6, and 10.5 meV. The first three 
peaks are globally consistent with the data from MARI 
[Fig. 3(b)]. However, the final peak, (lid), is very weak 
and does not seem to be present in the MARI data, so 
its existence must be said to be questionable. 



Figure 4(a) shows the inelastic spectra obtained on FO- 
CUS at temperatures of 1.5 K, 10 K, and 20 K, after 
integrating over a group of detectors such that the Q- 
range is 1 < Q < 3 A -1 . A vertical offset is applied 
between the curves for clarity. The results provide a sig- 
nificant improvement in the understanding of transitions 
(la) and (lb): these have energies of huJi a « 5.3 meV and 
hujji, w 6.2 meV at 1.5 K, consistent with the MARI data 
(taken at 8 K). As the temperature is increased there is 
a shift of both peaks towards lower energies, and some 
broad scattering develops below 5 meV. This is shown 
clearly by comparison of the data taken at 20 K and 
1.5 K (solid line), for which the same vertical offset is 
applied. The shift to lower energies arises from progres- 
sive population of the higher Ms sublevels in the S = 10 
ground-state multiplet, which results in a large number of 



E. Q-dependence 

We consider next the Q-dependence of transitions (la) 
and (lb). The integrated intensity is shown in Fig. 5 for 
the two transitions, with the FOCUS data represented 
by solid circles and squares, and the equivalent IN4 data 
by open symbols. Each integrated intensity point was 
obtained by fitting the peaks using a Gaussian function 
whose center, width, and intensity were all allowed to 
float. 

For both peaks, the intensity passes through a max- 
imum at Q w 1.45 A -1 , and decreases progressively at 
higher Q values. This Q-dcpcndcncc is typical of mag- 
netic scattering, where the intensity is expected to follow 
the magnetic form factor \F(Q)\ 2 , a rapidly decreasing 
function above Q ~ 2-3 A -1 , modulated by a structure 



7 




Energy Transfer (meV) 

FIG. 6: (a) INS spectra obtained on MARI at E % = 12 meV 
and T = 8 K for different values of Q. (b) INS spectra from 
IN4 at A = 2.2 A and T = 1.5 K for different values of Q. 

form factor I{Q) which depends on the exchange cou- 
plings and the exchange connectivity within the cluster. 50 
In the simplest exchange-coupled system, a dimer, the 
structure form factor is given by 50 

where R is the inter-atomic distance in the dimer. Using 
the dimer expression as an approximate indication, the 
maximum observed at Q w 1.45 A -1 corresponds to R w 
2.8 A, a value close to the shortest Mn-Mn separation 
in Mni2-acetate, d = 2.77 A. By comparison, the in- 
dependence of the ground-state transition at 1.25 meV 
is peaked at Q « 0.95 A" 1 , 6 consistent with its different 
physical origin. 

Figures 6(a) and (b) show the behavior of the INS 
spectra obtained respectively on MARI and IN4 over the 
full Q-range of each instrument. In both cases, the low- 
energy transitions (I) and (II) discussed above are lost 
rapidly in the incoherent phonon background as Q in- 
creases. The Q-dependence of peaks (III) and (IV) is 
more difficult to determine due to their weak intensity, 
but it is consistent with transitions of magnetic origin 
in that the intensity decreases with increasing Q. It is 



clear from both panels of Fig. 6 that the 27 meV peak is 
more pronounced at small Q values. A general increase 
of scattering is observed as Q is increased for all energies 
between 10 and 45 meV, which unfortunately impedes a 
complete analysis. 



F. Assignment of Magnetic Peaks 

Several magnetic peaks have been observed in four 
groups with energies faoi as shown in Table II. For transi- 
tions (la) and (lb), a detailed analysis may be performed 
of their dependence on Q and on temperature; for the 
higher transitions, only a more qualitative treatment is 
possible with the present data. 

To determine the nature of these transitions, we begin 
by noting that due to the negative axial zero-field split- 
ting only the Mg — 10 component of the S = 10 ground 
state is populated (over 99.5%) at 1.5 K, and therefore 
that only INS transitions from the state (10, ±10) can be 
observed. From the selection rules in Sec. IIIB, it follows 
that only states with 5 = 9, 10, and 11 can be excited by 
INS, a clear example of which is the observation 7 of the 
transition from £(10, ±10) to £(10, ±9) within the S = 
10 ground-state multiplet at fao « 1.25 meV [peak (GS1) 
in Fig. 3(c)]. This result leads to the accurate determi- 
nation of the axial anisotropy parameter for the S = 10 
ground state, D® = —0.0566 meV. 

The assignment of spin states for the excited energy 
levels depends both on experimental observation and on 
certain assumptions concerning the axial anisotropy pa- 
rameters for these states. We present first our treatment 
of anisotropy parameters, followed by the experimental 
results and a discussion of the consistency of this ap- 
proach, which leads to the spin assignments shown in 
Table II. As noted in Sec. II B and discussed in greater 
detail in Sec. IV, given the nature of the data and the 
energy scales involved we restrict our considerations to 
the level of anisotropy parameters of the total spin. The 
canonical energy-level splitting for a given S state may 
be expressed as 

E(S, M S ) = D S {M 2 S - ±S(S + 1)) , (3) 



TABLE II: Summary of energies and spin states of the ob- 
served magnetic transitions. 



1 


hu)i 


Properties 


la 


(5.3 ±0.1) meV 


S = 9, M s = ±9 


lb 


(6.2 ±0.1) meV 


S = 9, M s = ±9 


Ha 


(8.0 ±0.1) meV 


S = 9, M s = ±9 


lib 


(8.7 ±0.1) meV 


S = 9, M s = ±9 


lie 


(9.7 ±0.1) meV 


S = 9, M s = ±9 


lid 


(10.5 ±0.1) meV 


not confirmed 


III 


(27 ± 1) meV 


S = 9, 10, or 11 


IV 


(31 ± 1) meV 


may be magnetic 
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where D$ is the axial anisotropy parameter of the clus- 
ter in spin state 5. We will not attempt to compute the 
anisotropy parameters D$ from vector coupling of the in- 
dividual single-ion anisotropy terms, but for a qualitative 
interpretation of our observations will use related consid- 
erations to deduce their values relative to the measured 
ground-state axial anisotropy parameter D\ Q . In the ab- 
sence of experimental information on this point, we begin 
by assuming that D$ is the same for all excited states of 
the same spin 5, and return below to a more detailed 
consideration of this approximation. 

A qualitative picture of the relative values of D s for the 
lowest-energy states of each spin 5 in the Mni 2 -acetate 
cluster may be obtained from a crude model based on 
two observations. From our approximation to the crys- 
tallographic symmetry all single-ion anisotropy terms are 
identical for both types of ion, and from the Jahn-Teller 
distortions these are very much greater for the 8 Mn 3+ 
sites than for the 4 Mn 4+ sites. A vector-coupling cal- 
culation of the spin state anisotropy performed only for 
the Mn 3+ crown sites gives a systematic decrease of \D S \ 
with decreasing 5, and no change in the sign. In this 
elementary approach the anisotropy parameter deduced 
for 5 — 9 states, assumed to be composed of 3 Mn 3+ - 
Mn 3+ pairs with 5 = 4, one such pair with 5 = 3, and 
an antiparallel 5 — 6 core, is approximately 8% smaller 
than that for the 5 — 10 ground state. From a compar- 
ison with the coefficient D® Q = —0.0566 meV extracted 
from INS data in Ref. 7, this result is fully consistent 
with the value Dg w —0.049 meV estimated by fitting 
the magnetization to a similar dimerized model. 51 

For the uniaxial anisotropy parameter of the lowest 
excited 5 = 11 states, the simplest argument is to view 
these as a fully aligned S — 16 crown of Mn 3+ ions, with 
changes of spin state from the S = 10 ground state occur- 
ing on the core Mn 4+ ions, which have a very small single- 
ion anisotropy, whence w D® . The vector-coupling 
scheme above gives a value for somewhat larger than 
Di Q , but there is no qualitatively significant change in 
size or sign. The most important difference between 5 — 
9 and 5 = 11 excited states for our purposes is that the 
latter should show three excitation branches correspond- 
ing to the (S,M S ) final states (11, ±11), (11, ±10), and 
(11, ±9). These should be similar in intensity (differing 
by less than one order of magnitude) and should have a 
total energy span AE n = 40£>° 1 [see Eq. (3)]. 

We now demonstrate that the temperature-dependence 
of magnetic peaks (I) is consistent with transitions from 
the ground state (10, ±10) to excited states with 5 = 9 
and Ms — ±9. The Q-clependence of the magnetic tran- 
sitions at hu>i a — 5.3 meV and hu>n, = 6.2 meV is shown 
in Fig. 5. The decrease of intensity at low Q values is 
typical of transitions with AS = ±1, which have zero in- 
tensity at Q — 0. 50 In contrast, transitions with AS = 
arc strong at Q = 0. Further, as the temperature in- 
creases (Fig. 4) there is a growing scattering intensity 
below hu>i a . This suggests that additional magnetic tran- 
sitions are emerging because the increasing temperature 
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FIG. 7: (a) Detail of the low-energy levels inferred from our 
INS data (see text). After subtraction of the anisotropy shift, 
the lowest S = 9 states appear at 4.0 meV and 4.9 meV 
above the S — 10 ground state, (b) Energy levels observed 
by INS over the full range investigated. The spin values and 
anisotropy parameters of the higher levels are not known pre- 
cisely. 



causes thermal population of more Mg sublevels of the 
S = 10 ground-state multiplct. From the value and sign 
expected for the parameter Dg, such a situation is possi- 
ble only if the excited states corresponding to transitions 
(la) and (lb) are two separate S = 9 states. 

By applying Eq. (3), the energies of transitions from 
the ground state (10, ±10) to S = 9 spin states with 
Ms = ±9 can be written in the form 



hw a = £(9, ±9 
= E a + 



£(10, ±10) 



190 1 r>° I 
— l-^io I 



-51|D 9 °| 



(4) 



where E a is the energy of the 5 = 9 state in the absence 
of anisotropy. The multiplct structures are illustrated in 
Fig. 7(a), which shows also that the energy spread within 
the ground-state multiplet (lOO-D ^) is larger than that 
within the 5 = 9 manifolds (81 Dg). As temperature in- 
creases, the state (10, ±9) and other, higher states within 
the 5 = 10 multiplct are populated progressively and give 
rise to more INS transitions, which will necessarily occur 
at lower energies. 

The measured changes in intensity can be fitted with 
semi-quantitative accuracy using the value Dg = —0.045 
meV, and in fact our data suggest that the result of 
Ref. 51 should be regarded as an upper bound. With 
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these cluster anisotropy parameters, the energetic shift 
induced in the transitions between the lowest 5 = 10 
and 5 = 9 states is then 

fiw a «£ a + 1.29meV. (5) 

By contrast, if the excited states were of spin S = 11, 
for a reasonable value of one would expect the op- 
posite process (stronger scattering at higher energies), 
which is clearly precluded by the present data. Further, 
an explanation based on one or more 5 = 11 state(s) 
would require either a third excitation peak or that 
be very much smaller than D® so that all of the three 
peaks could be contained within the resolution width of 
the measurement; both of these possibilities may be ex- 
cluded. 

The situation is less clear regarding transitions (Ila-d) 
between 8 and 10.5 mcV: it is evident from their evolu- 
tion with temperature and with Q that these are mag- 
netic excitations. The intensity peaks undergo a small 
shift towards lower energies as the temperature is in- 
creased [Fig. 4(b)], while their wavevector-dependcncc 
(not shown) appears to have a weak maximum at in- 
termediate Q. In principle these transitions could corre- 
spond to final states with S = 9, 10 or 11. However, the 
combined T- and Q-dcpendence suggests that an 5 = 10 
state is unlikely. In distinguishing between 5 = 9 and 
11 states one may invoke both the weak T-dependence in 
combination with expected values for and the fact 
that the magnetization data of Refs. 38,39 do not appear 
to show a transition to an 5 = 11 state below a field of 
200 T, implying that the first candidate 5 = 11 peak in 
the INS data would be (III). Thus from the qualitative 
similarity in behavior of peaks (II) with peaks (I), we 
assign peaks (II) to transitions with 5 = 9 final states. 
In Sec. IV we will show that a numerical analysis of the 
Mni2-acetate exchange model confirms this assignment, 
and can be used to exclude the possibility of a set of 5 = 
11 transitions with significantly different from D® . 

We comment briefly on the assumption that the clus- 
ter anisotropy parameters D$ have the same value for 
all states of the same total spin 5. We have used a sim- 
ple model based on the classical spin arrangement of the 
5 = 10 ground state and the strong difference in single- 
ion anisotropies between Mn 3+ and Mn 4+ ions to offer a 
consistent interpretation of all of the measured INS data. 
A small extension of this scenario is to consider the clus- 
ter as 4 Mn 3+ -Mn 3+ -Mn 4+ trimer units, constructing the 
ground state from 4 5 = 5/2 trimers and excited states 
as perturbations about this. Under the same assumption 
concerning single-ion anisotropies, one obtains a range 
of values of Dg for different 5 = 9 states from 5% to 
100% of D° 1Q ; for excited 5 = 10 and 5 = 11 states, the 
values of D w and Dn range from 40% to 80% of D® . 
While this analysis is far from systematic, it illustrates 
three important qualitative features. First, it is difficult 
to find states with D > D° , a result which is not sur- 
prising given the fully aligned crown spins in this state, 
and helps to verify the assumptions underlying the in- 



terpretation. Second, the sign of D is always the same, 
suggesting that the temperature-dependence of the ob- 
served peaks can be used as a reliable indicator of the 
magnitude of Dg- Finally, this magnitude may in fact 
depend rather strongly on 5; however, the constant sign 
is sufficient to ensure that the conclusions we have drawn 
from our INS data are robust against such quantitative 
changes. 

In conclusion, the Q-dependence of transitions (la) and 
(lb) corresponds to transitions from the 5 = 10 ground 
state to 5 = 9 final states, and not to 5 = 11 or other 
5 = 10 states. While a similarly quantitative analysis 
is not possible for the transitions (Ila-d), we have ob- 
tained sufficient evidence to conclude that these are also 
of 5 = 9 character. For transitions (III) and (IV) a de- 
tailed investigation is excluded, and, on the basis of the 
general intensity decrease with increasing Q, we state 
only that these are of magnetic origin. A summary of 
the magnetic excitations identified by INS, their ener- 
gies, and their 5 values is presented in Table II, and a 
graphical representation is provided in Fig. 7(b). 

IV. MODEL ANALYSIS 

The microscopic Hamiltonian for exchange interactions 
between individual Mn ions in the Mni2-acetate cluster 
(Fig. 2) may be expressed as 

4 

H = J\ S$j-2 ■ 5*3j-i + J3 S$j-2 ■ S^j + i 

J=l 

+ (^2 Ssj-2 + J4 53j_l) • (53j_3 + S^j)} 

+J 3 (S 1 -S 7 + S 4 -S 10 ), (6) 

with periodic boundary conditions i = i + 12. Sites i = I, 
4, 7, and 10 represent the core 5 = 3/2 spins on the Mn 4+ 
ions, with mutual interactions J3, and the other sites cor- 
respond to the 5 = 2 crown spins (Mn 3+ ), which have 
exchange interactions J\ and J2 with different core spins 
and J4 with each other around the outer ring (Fig. 2). 
As stated in Sec. II A, we neglect other possible exchange 
paths, and assume despite the small structural distor- 
tions that the system is close to full fourfold rotation and 
reflection symmetry, and thus that J 2a = -hb = J2 and 
Jia = Jib = Ji- We will return later to a quantitative 
statement on this latter point. For the numerical anal- 
ysis to follow we will not consider single-ion anisotropy 
terms around individual ions, which are instead included 
through the term in Eq. (3) for the spin states of the clus- 
ter. We will also discard any of the possible higher-order 
anisotropy terms, such as those of the form 5 4 and 
\B\ (5 4 + 5i). These approximations may be justified 
on the grounds of the small energy scale of the terms 
involved in comparison with the scale of the magnetic 
excitation spectrum (Sec. III). 

The aim of this section is to determine a set of ex- 
change constants { Ji} which is capable of explaining the 
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magnetic excitation spectrum measured in Sec. III. Spec- 
tra of low-lying magnetic excitations have been tabulated 
by Raghu et a/. 43 by considering a variety of parameter 
sets. It is evident from this analysis that none of the 
sets considered yields sufficiently many excitations with 
appropriate total spin below an energy of approximately 
300 K (30 meV) to explain the INS results. Similarly 
the excitation spectra computed for the parameter set 
of Ref. 44 also fails to provide a sufficient number of 
low-lying excitations to account for the new INS data 
(although we note here that these contain rather more 
levels in the relevant energy range than previous authors 
were aware of). While the parameter set of Ref. 46 is 
not inconsistent with the excitation spectrum, we will 
find in Sec. IV B that it cannot account for the magnetic 
susceptibility data. Hence the parameter sets discussed 
in Sec. II D are all incompatible with the experimental 
data. We therefore begin with an unbiased determina- 
tion of an appropriate set of exchange constants {Ji} by 
considering the ground state of the system using the full 
available parameter space. Candidate parameter sets are 
tested by computing the magnetic susceptibility, which 
we compare with our own measurements in Sec. IV B. In 
Sec. IV C we comment briely on the high-field magnetiza- 
tion data. We return to the issue of magnetic excitations 
in Sees. IV D and IV E, where we present a theoretical 
analysis of the energy spectrum of the cluster. 

A. Numerical Methods 

The numerical results are obtained by two methods. 
We have performed a systematic high-temperature series 
expansion to compute the magnetic susceptibility 

T X (T) = C + C x + C 2 2 + C 3 3 + C 4 /3 4 + . . . , (7) 

where = l/k^T with /cb the Boltzmann constant. For 
illustration, the first four coefficients C n for the cluster 
geometry and ionic spins of Mni 2 -acetate are given by 

C = 21, 

d = -20Ji - 40J 2 - f J 3 - 64J 4 , 
C 2 = -5J? + 80JiJ2 + 55Jf + 75JiJ 3 + 150J 2 J 3 
+ if J 3 + 80J i J 4 + 160J 2 ^4 + H2J|, 

c 3 = ^-J\ + 20J 2 J2 — soji j 2 + ^fji ~ ^T'JiJz 

-300JiJ 2 J 3 - 325J 2 2 J 3 - £|£ Ji J 3 2 - £|£ J2J3 2 

- W J 3 + 20J i j4 _ 380Jl J 2 Ji - 280J 2 2 Ji 
-300Ji J 3 Ji - 6OOJ2 J 3 Ji - 140Ji Jj 
-280J 2 J 4 2 - W J l (8) 

One observes that Co is a constant, C\ a weighted sum 
of the cluster interactions, C2 a sum of certain combina- 
tions of squares, and so on. We have continued the series 
to 8th order, where the expressions are best handled by 
symbolic computation methods, 52 and have included the 



full Mni2 cluster rather than employing a conventional 
linkcd-cluster expansion. 

The second method we employ is the calculation of 
low-lying excitations for a given set of input interaction 
parameters {Ji} by exact diagonalization (ED) of the 
cluster Hamiltonian (6) using the Lanczos procedure. We 
have used the conservation of S z , the z-component of the 
total spin, as well as the spatial symmetries of the cluster. 
Conservation of S z is particularly important in reducing 
the dimension of the Hamiltonian matrix and is one of the 
reasons for which we do not include single-ion anisotropy 
terms in Eq. (6). Spatial symmetries of the states are 
denoted by a "momentum" k such that a state acquires 
a factor exp(ifc) under a 90° rotation of the model for the 
Mni2 cluster shown in Fig. 2. We have used in addition 
the reflection symmetries of the cluster for k = and 
k = it. For the description of the INS results, and to 
test whether the spin of the ground state is S = 10, one 
must consider at least the sectors with S z > 9. The 
largest matrix dimension is then 324 908, which occurs 
for S z = 9, k = it/2. This dimension is sufficiently small 
that rapid diagonalization by the Lanczos procedure is 
possible on a modern personal computer. Indeed, fast 
diagonalization is a necessary condition for an analysis 
requiring the consideration of many sets of parameters 
{ Ji}. Finally, using our optimal parameter set we have 
performed further computations with S z = to ensure 
that the spin of the ground state is indeed 5* = 10. For 
S z = the largest dimension to be considered is 1 073 763 
when using spin inversion symmetry. 

The condition of a ground state with spin S = 10 al- 
ready sets a strong constraint on the ratios of the ex- 
change interactions, and, for example, may be used to 
exclude the parameter sets proposed in Refs. 16 and 
45. We have computed low-energy spectra for approx- 
imately 1500 independent parameter sets which satisfy 
the S = 10 ground-state condition. Many of these sets 
were considered in Ref. 44, but we have performed inde- 
pendent computations for all parameter sets and added 
further points in regions indicated by the susceptibility 
calculation as potentially relevant. We note that the 
S = 10 condition is invariant under rescaling of the in- 
teraction parameters, whence one obtains a line in the 
four-dimensional parameter space from each independent 
diagonalization. 



B. Magnetic Susceptibility 

The behavior of the static susceptibility of a magnetic 
material at high temperatures may be used in a system- 
atic manner to extract information concerning the ex- 
change interactions within the system (see for example 
Ref. 53). The simplest example of this process is that 
the first correction to paramagnetic behavior in a FM or 
AFM system gives the Curie- Weiss temperature, which is 
directly proportional to the sum of the couplings to each 
individual spin. In Mni2-acetate the application of this 
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FIG. 8: Static magnetic susceptibility measured with an ap- 
plied field of 0.1 T (points) and the best fit by obtained by 
high-temperature series expansion (solid lines). The param- 
eter set used is given in eq. (10). The solid lines are ob- 
tained from different approximations to the 8th-order high- 
temperature series, and represent the error estimate for the 
extrapolation procedure. 

analytical approach is complicated by two factors: first, 
the susceptibility is still far from its high-temperature 
limit when the sample decomposes above 300 K, and 
second, a number of the coupling coefficients to be de- 
termined may be of similar magnitude. 

We begin by showing our measurements of the mag- 
netic susceptibility of Mni 2 -acetate, which are given by 
the symbols in Fig. 8. This measurement was performed 
on a 50 mg powder of the fully deuterated sample, over a 
temperature range of 1.8-300 K and with an applied field 
of 0.1 T, using a Quantum MPMS XL-5 magnetometer, 
and the data were corrected for diamagnetic contribu- 
tions. Our data is very similar to that obtained in Ref. 54, 
with the exception of the dip at the lowest temperatures, 
which is significantly sharper in our 0.1 T measurement 
because the spread of sublcvcls in the S = 10 ground 
state is narrower than in an applied field of 1 T. 54 We 
comment that the presence of a small component of a 
faster-relaxing species of Mni 2 -acetate molecules, esti- 
mated at 3-4% in Ref. 7, would be visible only at the 
lowest temperatures. 

The measured susceptibility is then fitted in the form 

oo 

T X {T) = Y J C n T- n . (9) 

n=0 

Making use of the g- values obtained for Mni 2 -acetate by 
EPR measurements on a powder sample, 21 we begin by 
fixing the value of the leading coefficient as Co = 30.5 
emu K mol -1 . The higher coefficients C ni n > 1, are 
then determined from a fit to the high-temperature part 
of the data, from which one may obtain an estimate of 
the uncertainty in the coefficients from their variation 
with the temperature range used for the fit. 

The interaction parameters J\, J 2 , Js, and J4 are then 
fitted by comparing the coefficients C n and C n in Eqs. (7) 
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FIG. 9: Suitability of parameter sets giving a ground state 
of spin S — 10 for fitting the high-temperature susceptibility 
data: (a) as a function of Ji and J2 for all values of J3 and 
J4; (b) as a function of J3 and J a for all values of Ji and J2. 
Darker shades indicate superior agreement with x(T). 



and (9). Because this is an overconstrained problem, 
and in order to satisfy the condition that the ground 
state be S = 10, we have optimized a suitability function. 
This function was in turn defined as a sum of differences 
of appropriately rescaled coefficients C n and C n , 1 < 
n < 6, normalized by the estimated uncertainties in the 
latter coefficients. The single most important piece of 
information which can then be extracted is C\, which is 
given by a weighted sum of the interactions (see Eq. (8)) 
and hence sets the overall energy scale of the problem. 

The result of this analysis for fixed Cq = 30.5 emu K 
mol -1 is shown in Fig. 9. Each point in this figure is a 
projection of a point in the four-dimensional parameter 
space yielding an S = 10 ground state. Fig. 9(a) shows a 
projection of the four-dimensional parameter space onto 
the (Ji, J 2 ) plane and Fig. 9(b) onto the (J3, J4) plane. 
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The suitability of the parameter sets for describing the 
measured susceptibility xC^O is shown by shading of the 
points, where darker shades denote higher and lighter 
shades lower levels of agreement. In Fig. 9(a) we observe 
two maxima of suitability. The left maximum gives rise to 
5 = 11 excitations with energies of approximately 250 K 
(« 21.5 mcV), an energy range in which from Sec. IIIF 
no such excitations are expected, and thus we focus on 
the other maximum. We conclude that the requirements 
for an S = 10 ground state, for matching the susceptibil- 
ity x(T), and for S = 11 excitations of sufficiently high 
energies, constrain the exchange parameters (to within 
uncertainties of order 3-4 K), to the values J\ w J2 w 61 
K, J 3 w 7.5 K, and J 4 « 10 K (Fig. 9). 

In this region of parameter space, we have further op- 
timized a definitive set of exchange interactions by con- 
sidering the magnetic excitations, which are discussed in 
detail in Sec. IV D. The primary constraint set by the 
INS measurements presented in Sec. Ill is the require- 
ment that the lowest S = 11 excitation have an energy 
of at least 285 K (w 25 meV). This condition moves the 
optimal parameter set away from the minima of the suit- 
ability function found in Fig. 9 for the illustrative fitting 
procedure with fixed Co, yielding the set 

Ji = 67.2 K (5.79 meV), J 2 = 61.8 K (5.33 meV), 
J 3 = 7.8 K (0.67 mcV), J 4 = 5.6 K (0.48 mcV).(10) 

We will focus henceforth on this parameter set, and be- 
gin by justifying it as our final choice. For this purpose 
we present a direct comparison of the magnetic suscepti- 
bility obtained with the parameters of Eq. (10) and the 
experimental data. Instead of fixing the precise value 
of the zeroth-order coefficient we use this as a further 
fitting parameter, obtaining the slightly modified coeffi- 
cient Co = 29.5 emu K mol -1 . This value corresponds 
to an effective average (/-factor g e g = 1.935, and the 
modification can be interpreted as arising from single-ion 
anisotropies which are not otherwise present in the anal- 
ysis (Sec. IIIF). In addition, we now use all eight orders 
of the high-temperature expansion, two more than were 
used in the calculation of the suitability function shown 
in Fig. 9. The lines in Fig. 8 show several Pade approxi- 
mants to the series obtained from the optimal parameter 
set. Without entering into the technical details of this 
procedure, we note only that the different Pade approxi- 
mants to the series allow an estimate of the uncertainty 
in the extrapolation of the high-temperature expansion. 
The divergence of the different approximants below 150 
K, and the departure from the data at low temperatures, 
are to be expected in a high-temperature series expan- 
sion. Within the uncertainty of the extrapolation, we 
find good agreement with the data over the entire tem- 
perature range from approximately 80 K to 300 K. In the 
high-temperature regime which is relevant for our anal- 
ysis, the effect of a small concentration of faster-relaxing 
species, i.e. of molecules with slightly different exchange 
constants, would be negligible both in the suitability of 
the fit and in comparison with the other sources of un- 
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FIG. 10: Calculated magnetization curve for a Mni2-acetate 
cluster with exchange parameters given in eq. (10). 



certainty We stress that the overall agreement is in fact 
better than that provided by the maxima in Fig. 9 un- 
der the constraint Co = 30.5 emu K mol -1 . We conclude 
that the optimized parameters of Eq. (10) provide a good 
description of the experimental data for the magnetic sus- 
ceptibility. 



C. High-field Magnetization 

We comment here that the constraint on the overall 
magnitude of the interaction parameters provided by fit- 
ting the susceptibility requires a new interpretation of 
the high-field magnetization data. 38 ' 39 The fixed value of 
the weighted sum of exchange constants C\ has a direct 
correspondence to the saturation field of the system, at 
which all bonds must be polarized ferromagnetically. The 
predicted T — magnetization curve, by which is meant 
here the spin 5 ~ M/g/iB as a function of applied field 
H, is shown in Fig. 10 for the parameters of Eq. (10), 
using the same value <7 e ff = 1-935 as for the susceptibility 
to express the magnetic field H in units of T. The low- 
field data is consistent with that of Refs. 38,39, in that 
the S = 10 ground state persists up to fields in excess 
of 200 T. The relatively large number of narrow steps in 
the calculated magnetization may then also be compared 
with the regularly spaced steps in the data. The location 
of the first step, corresponding to the transition between 
S = 10 and S = 11 ground states, may not be under- 
estimated from our computation by more than 10-20 T. 
The optimal susceptibility fit therefore implies peaks in 
the dM/dH data at fields smaller than the lowest peaks 
assigned and used to determine the exchange parameters 
{Ji} in Ref. 44. However, we believe that the raw exper- 
imental data 38 ' 39 do not in fact exclude further peaks in 
the region between 200 and 300 T. 

The uncertainty due to errors in the exchange con- 
stants may accumulate to a value on the order of 100 
T for the saturation field, which is the transition be- 
tween S = 21 and 22 ground states, but larger devia- 



13 



tions become increasingly unlikely given the constraint 
on the overall energy scale obtained from the susceptibil- 
ity. Hence the single largest peak in the dM/dH data of 
Refs. 38,39, which appears at 520 T, occurs at or above 
the predicted saturation. Indeed, none of the magneti- 
zation curves which we have computed yield indications 
for a single step which would be substantially more pro- 
nounced than all of the other steps. Given the reliabil- 
ity of a susceptibility measurement, and the destructive 
nature of the explosive compression technique involved 
in the high-field magnetization measurement, we suggest 
that the former is more representative. In this interpreta- 
tion, all signals at and beyond the largest peak at 520 T 
in the dM / dH data would correspond to times when the 
sample or other components of the experimental appara- 
tus are already disintegrating. 



D. Excitation Spectrum 

We turn next to an analysis of the excited states. Ta- 
ble III shows the lowest nondegenerate magnetic energy 
levels for a cluster with the optimized set of exchange in- 
teractions (10). These levels are labeled by their spin and 
spatial symmetry sectors. Qualitatively, the parameter 
values can be seen to provide a number of low-lying S = 
9 states, as required for comparison with INS data [tran- 
sitions (I) and (II)], without permitting the existence of 
an S — 11 state or additional S = 10 states below a 
significantly higher energy, whose value is in agreement 
with the limits set by INS [transition (III)] and high-field 
magnetization measurements. The fact that five nondc- 



TABLE III: Magnetic energy level spectrum obtained for the 
first 17 nondegenerate states with S > 9 of a Mni2-acetate 
exchange model with the parameters of Eq. (10), classified 
according to spin state, level degeneracy, and spatial sym- 
metry sector. No energetic correction is applied for uniaxial 
anisotropy of the cluster. 



Spin S 


Energy [K] 


Degeneracy 


Symmetry k 


10 


0.00 


1 
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28.48 


2 


±7r/2 
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44.47 
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7T 
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91.46 


1 
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119.67 
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±7r/2 


9 


159.61 
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7T 


11 


285.58 
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304.45 


1 





9 


304.65 


1 


7T 


9 


306.73 


2 


±7r/2 


9 


307.14 


1 


7T 


10 


311.11 
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FIG. 11: Comparison with INS data of the location of cal- 
culated magnetic excitations, obtained using the optimal pa- 
rameter set (10). An anisotropy correction is applied by us- 
ing Eq. (3) with a single D value for all (ground and excited) 
states of each total spin 5 = 9, 10, or 11. The experimen- 
tal data was obtained on MARI at (a) low energy transfer 
[Figs. 3(b) and (c)] and (b) high energy transfer [Fig. 3(a)]. 



generate S = 9 states are found in the low-energy man- 
ifold may be justified by straightforward considerations 
based on a spin- wave description (Sec. IV E), and can be 
used as an aid to experimental interpretation (Sec. IIIF). 
By considering the symmetry of the S = 9 states one may 
also observe which of these would be sensitive to a break- 
ing of fourfold cluster symmetry, and thus by compari- 
son with the width of the observed INS peaks establish 
an approximate upper bound on the extent of any such 
departure from symmetry (Sec. IV E). 

A quantitative comparison between the numerical re- 
sults and the positions of the excited levels observed 
by INS requires that the single-ion anisotropy terms be 
taken into account. As stated above, we have included 
these terms only at the level of an effective uniaxial 
anisotropy acting on the cluster spin [Eq. (3)]. The re- 
sults of this exercise are shown in Fig. 11. With regard 
to position, type, and number of groups of peaks, the 
properties of the model are in good accord with the mea- 
sured data sets. Levels (la), (lb) and (Ha) deviate by no 
more than 20 K (« 2 meV) from the theoretical results 
while the discrepancy for levels (lib) and (lie) is some- 
what larger. In this regard, we note that the anisotropy 
shift is on the order of 15 K (k 1.3 meV) and that rather 
small changes (of order 1 K) in the values of any individ- 
ual coupling parameter Ji may cause large (of order 10 
K), albeit not entirely uncorrelated, changes in the posi- 
tions of the different excited states. Further optimization 
would require a full analysis of the uniaxial anisotropy at 
the single-ion level. Although computation of spectra for 
a number of selected parameter sets with the inclusion 
of single-ion terms is possible with current computing 
technology, it is a demanding numerical task due to the 
fact that S z is not a conserved quantum number in this 
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case. As a consequence the full optimization of param- 
eters would no longer be possible, and for this reason 
we do not pursue single-ion terms further here; the en- 
ergy scale of the resulting mixing and splitting of levels 
in different manifolds remains small on the scale of the 
variability present in the exchange interactions. We thus 
focus on the robust features which may be extracted at 
the current level of refinement, namely a reproduction 
of the properties of the peak groups (I), (II), and (III) 
(Fig. 11). 

In addition to the position of the magnetic excita- 
tions it is important also to compute their intensities 
to ensure that transitions to all of the calculated lev- 
els are allowed. We have computed the matrix elements 
(g(10, 9)\Sf\e n (9, 9)), which by spin rotation symmetry 
are related directly to the elements of the operator S~ 
appropriate for the spin-flipping action on the 5 = 10 
ground state \g(10, 10)) of a neutron-scattering process, 
for all sites i = 1, . . . , 12 of the cluster. These elements 
may be used in a calculation including the structure fac- 
tor of the molecule to obtain the transition matrix ele- 
ment at any wave vector Q. In Table IV we show only 
the incoherent sum (no phase factor) of the squares of the 
site matrix elements: this is sufficient to show that none 
of the levels is excluded for symmetry reasons, and that 
all may be expected to have similar weights. However, 
as a result of the degeneracy of the pairs (|ei), \e2)) and 
(1 65), for INS purposes the relevant quantity would 
be the sum of these two weights. One may then expect 
two peaks whose integrated intensity is double that of 
the other three, a result certainly not inconsistent with 
the INS data. 

We conclude this section by computing the spin distri- 
bution of the ground-state wave function. Table V shows 
the values of (5?) 2 on the core 5 = 3/2 spins (i = 1, 4, 
7, 10; Fig. 2), on the 5 = 2 crown spins coupled to one 
core spin by exchange parameter J\ (i — 2, 5, 8, 11), and 
on the 5 = 2 crown spins coupled to two core spins by 
exchange parameter J 2 (i = 3, 6, 9, 12). The remark- 
able feature of these results is that approximately 90% 
of the spin weight of the crown spins is in the classical 
ferrimagnetic state invoked in Sec. II B to justify the S 
— 10 ground state. This value drops only to 80% for 
the core spins, whose mutual coupling is weak but frus- 
trating. The degree of overlap with the classical state 



is confirmed by considering the S z correlations between 
neighboring incquivalcnt sites (Table V). 

The classical alignment is ensured by two factors. The 
first is the dominance of the exchange couplings J\ and 
J2, which forces all of the core spins be antiparallel to 
the crown spins, and therefore mutually FM. The weaker, 
frustrating J3 and J4 bonds then cause only minor devi- 
ations from the classical state. The second is the near- 
equivalence of the interactions J\ and J2, to which the 
sensitivity of the system is shown by the size of the 
black region in Fig. 9(a). The effect of a departure from 
this equivalence is illustrated by the example of the J\- 
dominated parameter set originally proposed in Ref. 16, 
where sites i = 1, 2, 4, 5, 7, 8, 10, and 11 would be cou- 
pled to form effective 5=1/2 units, and only very small 
ordered spin components would be observed on these. 
Precisely because corrections of this form are small in 
Mni 2 -acetate, one may expect that a spin-wave descrip- 
tion of the excited states 55 is in fact meaningful despite 
the small size of the system. Finally, the relative lack of 
quantum mechanical fluctuation effects may also be as- 
cribed in part to the "large" values of the ionic spins (5 
= 3/2 or 2, as opposed to 5 = 1/2). 

E. Spin-wave Analysis 

We have found that the optimal parameter set for 
Mni2-acetate yields a ground state whose spin distribu- 
tion can be said to be "80-90% classical". This result, 
which may seem surprising for a relatively small, quan- 
tum system, is a consequence of the fact that the two 
dominant exchange interactions, J\ and J2, are unfrus- 
trated. The observation that an 5 = 10 ground state may 
be found only in a very small region of four-dimensional 
parameter space around this rather special limiting case 
is evidence that any significant frustrating interactions 
J3 or J4, or deviations from J\ « J2, would in fact de- 
stroy the classical ground state. The gap to the lowest 
excited states of the exchange model without anisotropy, 
measured as being approximately 4 meV (40K) in Mni 2 - 
acetate [Eq. (5), Fig. 7], may also be considered as a sen- 
sitive indicator of the proximity of any similar parameter 
set to a ground state of a different spin 5 ^ 10. 

One consequence of this quasi-classical nature is that 



TABLE IV: Incoherent sum of squared matrix elements M n = 
J2i I (5(10, 9)|Sf |e„(9,9)}| 2 between ground and excited states 
for the 7 low-lying S = 9 levels, obtained with the parameters 
of Eq. (10). 



TABLE V: Spin distribution of the ground state of Mni2- 
acetate, characterized by the values of (Si) 2 on all sites and 
by spin correlations between neighboring inequivalent sites, 
obtained with the parameters of Eq. (10). 
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Energy [K] 






0/(io, io)|(sy) 2 | 5 (io,io)) 


1.77542 


i = 1, 4, 7, 10 


1,2 


28.48 


0.25506 




3.63387 


i = 2,5, 8, 11 


3 


44.47 


0.21426 






3.55449 


i = 3, 6, 9, 12 


4 


91.46 


0.19241 






-2.40483 


i = 1, j = 2 


5,6 


119.67 


0.19485 




0/(io, io)|Sfs;| ff (io, 10)) 


-2.32373 


i = l, j = 3 
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159.61 


0.17911 
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i = 2, j = 3 



15 



300 p 3 § m 

A 



* 



200 





o 


ED, S=9 


* 


V 


ED, S=10 




A 


ED, S=1 1 




* 


spin-wave 



LU * 

O 




71/2 



k 

FIG. 12: Comparison between magnetic excitations calcu- 
lated by exact diagonalization (ED) and within a spin-wave 
approximation using the parameter set of Eq. (10). No uni- 
axial anisotropy is included. Higher S = 9 and 10 excitations 
not included within the 3 spin-wave branches enter only at 
the highest energy of the 5=11 branch. 



the magnetic excitations may be considered by conven- 
tional spin- wave theory. Within this treatment, first pre- 
sented in Ref. 55, the cluster is considered as a ring of 
4 unit cells each with 3 inequivalcnt sites. Two of the 
spin-wave branches correspond to reduction of the spin 
component S z = 2 on the crown sites, leading to an S 
= 9 state of the classical cluster, while the third corre- 
sponds to an increase of the component S z — —3/2 on 
the core sites, leading to an S — 11 state. The 4-cell ring 
has 3 inequivalent k points, where k is an effective wave 
vector around the ring, whence one expects 9 levels in 3 
branches: the ground state, 5 S — 9 states, and 3 S = 11 
states. The exact energy levels for any parameter set {Jj} 
may be computed by exact diagonalization, and in Fig. 12 
we show the comparison between these levels and the 
spin-wave approximation for the optimal Mni2-acetate 
set (10). The agreement is remarkably good for such a 
crude approximation to a small quantum system, again 
because of the special, "classical" set of exchange inter- 
actions found in Mni 2 -acetate. The authors of Ref. 55 
applied their approximate treatment to the parameter 
set of Sessoli et al., 16 for which the ground state deviates 
strongly from the classical paradigm, and unsurprisingly 
a correspondence between their excited states and the ex- 
act results is difficult to establish in this case. Returning 
to the interpretation of the INS data, the peak groups I 
and II (see Table II) can be identified respectively with 
the first and second branches of spin- wave excitations, 
all of which have spin 5 = 9 and provide the appropriate 
symmetry sectors (k in Table III). 

The success of the spin-wave approximation in the 
physical parameter regime may be used for two further 
purposes. One is to investigate the sensitivity of the com- 
puted excitations to changes in individual interaction pa- 
rameters without the need for a lengthy numerical cal- 
culation. We have performed this exercise in order to 
obtain the optimized parameter set (10), and also to ob- 



serve the relative changes of different excited levels with 
the parameters. Here the fact that certain excitation en- 
ergies change together suggests that a full, independent 
tuning of the levels to obtain exact agreement with the 
measured INS peaks may not be possible within a four- 
parameter exchange model, and that this would proba- 
bly require the reintroduction of microscopic anisotropy 
terms neglected in the current analysis. In this connec- 
tion we note only that single-ion anisotropy terms are 
expected to make the spin- wave branches less dispersive, 
which would improve the agreement with the INS results. 
The second purpose is to investigate the consequences of 
a breaking of fourfold cluster symmetry, if in fact one 
were present. In this case J 2a ^ J 2b and J^ a ^ J^, and 
the degeneracy of the k — ±ir/2 states in Table III would 
be lifted. With reference to Table IV, these are the de- 
generate level pairs (|ei), |e 2 }) and (|e 5 ), |e 6 }). From the 
fact that the INS data may already be fitted by 5 Gaus- 
sian peaks of rather similar widths [r = 0.77 meV (8.9 
K) in (I) and (II)] not far from the resolution limit, the 
spin-wave description can be used to set upper bounds on 
AJ 2 = \J2a — J2b\ and AJ4 = \ J^ a — Jib\ of approximately 
0.1 meV (IK). 



V. SUMMARY 

We have performed inelastic neutron scattering mea- 
surements up to high energies to identify and characterize 
the magnetic excitations of Mni 2 -acetate. We find that 
all of the lowest energy levels, which occur in two groups 
at 5-6.5 meV (60-75 K) and 8-10.5 meV (95-120 K), 
appear to be S = 9 states. There are approximately five 
such levels, and their intensities do not vary by more than 
a factor of 2-3. Higher levels, including the lowest 5=11 
states, are found at a significantly higher energy, namely 
the newly identified transition peak (III) at 27 meV (310 
K) , which is qualitatively consistent with the location of 
the first magnetization step. Taken together with the fact 
that the ground state has spin S = 10, the qualitative 
features of the measurements are already sufficient to re- 
strict the interaction parameters of the cluster almost to 
a unique set. All of the available experimental informa- 
tion is reproduced by a parameter set with J\ ~ J 2 ~ 5.5 
meV (65 K), while J3 and J4 are smaller than 1 meV (10 
K); refinement of the more robust quantitative features 
leads to the parameter set given in Eq. (10). 

Essential supplementary fitting information is pro- 
vided by the magnetic susceptibility, which we have com- 
puted by a systematic high-temperature series expansion. 
A correct reproduction of the high-temperature limit sets 
a constraint on the four interaction parameters which 
determines the overall energy scale of the couplings in- 
volved. In this context we note that the leading term 
may be expected from the cluster geometry and dominant 
couplings to be determined by the combination Ji + 2 J 2 , 
which is indeed very close to the value J\ = 216 K ex- 
tracted in Ref. 16 on the assumption of one dominant 
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interaction. In fact this constraint is difficult to reconcile 
with the high-field magnetization data, whose saturation 
field should be given by a similar combination of terms. 
On the basis of our measurements and calculations, the 
largest feature in the dM/dH curve of Refs. 38 and 39 
should correspond qualitatively to the saturation field. 

In the consistent parameter set all four of the ex- 
change interactions are AF. From qualitative consider- 
ations based on exchange and supcrcxchange processes 
these results are fully plausible, despite the proximity of 
certain y(i-oxo bridging angles to 90°. The quantitative 
similarity of J\ and J 2 in spite of the very different types 
of exchange pathway involved, as well as the small but 
positive value of J3, emphasize the difficulties inherent 
in performing ab initio calculations sufficiently accurate 
to reproduce the physical properties of complex magnetic 
systems. Even small deviations from these values are suf- 
ficient to change the spin of the ground state (Sec. IV B), 
to bring excluded higher-spin states into the low-energy 
manifold, or to return a high-temperature susceptibility 
well outside the limits imposed by the experimental mea- 
surement. 

Throughout our analysis we have focused on the quali- 
tative features which a candidate parameter set must re- 
produce. At the quantitative level there are two sources 
of error: the first is the uncertainty in the multiparame- 
ter problem of fitting the set {J,} to disparate (and occa- 
sionally conflicting) pieces of data which are themselves 
subject to measurement errors; the second is the terms in 
the Hamiltonian which we have neglected, primarily ex- 
change asymmetries, single-ion anisotropies and higher- 
spin interactions. Because of the first, it is manifestly 
not meaningful to ascribe parameter values to the second, 
and yet these would be required for a perfect quantitative 
account of an ideal magnetic excitation spectrum. Thus 
we have restricted our considerations to robust results 
and qualitative fitting, adequate at the level of groups of 
INS peaks. These conditions are in fact sufficient to iden- 
tify an unambiguous set of exchange constants. Although 
for reasons of simplicity we have focused our presention 
on the precise values given in Eq. (10), we have performed 



similar analyses for small modifications of this set. From 
the variability of exchange parameters between sets yield- 
ing a description of similar quality, we estimate an error 
bar of 2-3 K for each individual exchange parameter Jj . 

In conclusion, measurements of its magnetic proper- 
ties allow us to establish a definitive set of intramolec- 
ular exchange interactions for Mni 2 -acetate. The only 
consistent parameter set is 

Ji = 67±3K (5.8±0.3meV), 

J 2 = 62±3K (5.3±0.3mcV), 

J 3 = 8±3K (0.7±0.3meV), 

J 4 = 6±3K (0.5±0.3mcV), 

where the error bars denote the order of uncertainty for 
each Ji. These parameters provide both a straightfor- 
ward understanding of the ground state and low-energy 
spectrum, and a resolution of the conflicts in the existing 
literature. 
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